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ABSTRACT 


A  parabolized  Navier-Stokes  Code  (HYTAC)  is  used  to  pre¬ 
dict  the  supersonic,  laminar  or  turbulent  viscous  flow  about  arbi¬ 
trary  geometries  at  moderate  angles  of  attack,  and  the  results  are 
compared  with  those  from  other  viscous  and  inviscid  codes.  The  test 
cases  presented  in  this  paper  include  the  external  flow  over  a  su¬ 
personic,  blunt  nosed  inlet,  a  slab  delta  wing  and  an  elliptic  body 
at  various  Mach  numbers.  The  results  for  the  inlet  are  compared 
with  the  results  from  a  viscous  shock-layer  code  (VSL3D)  and  they 
show  good  agreement.  The  slab  delta  wing  surface  pressure  compares 
well  with  the  results  from  an  inviscid  code  (STEIN)  and  with  experi¬ 
mental  results.  The  surface  heat  transfer  for  this  case  is  in  good 
agreement  with  the  results  from  a  boundary  layer  code  using  stream¬ 
line  tracing  (TRACE)  and  with  experimental  results./, 


INTRODUCTION 


I” 

k 


The  design  and  successful  flight  of  high  speed  inlets  and  complex  lift¬ 
ing  and  reentry  vehicles  can  be  aided  by  the  use  of  accurate  and  efficient 
computational  fluid  dynamics.  In  this  paper,  results  of  computer  analyses 
of  three  bodies  of  complex  shapes  are  given.  The  first  body  considered  is  a 
supersonic  blunt  nosed  inlet,  designed  to  operate  at  Mach  numbers  from  4  to  7 
and  at  angles  of  attack  up  to  10  degrees  (Figure  1).  The  inlet  is  a  sphere- 
cone  with  an  afterbody  of  smooth  curvature.  The  second  body  is  a  slab  delta 
wing  with  80  degree  sweep  angle  (Figure  2).  This  wing  has  a  nose  radius  of 
0.5  inch  and  is  designed  to  operate  at  a  Mach  number  of  9.6  at  various  angles 
of  attack.  The  third  geometry  is  an  elliptic  body,  which  has  a  spherical 
nose  of  0.3  feet  radius  and  an  afterbody  with  elliptical  cross  section 
(Figure  3).  The  change  from  spherical  to  elliptical  cross  section  is  gradual, 
and  the  final  cross  section  of  2:3  ellipse  is  achieved  within  8  to  10  nose 
radii.  This  body  is  designed  to  operate  at  Mach  numbers  of  about  20.  Flight 
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performance  under  all  the  above  conditions  requires  a  thorough  understanding 
of  the  aerodynamic  forces  and  the  surface  heat  transfer  and  of  their  in¬ 
fluence  on  the  vehicle  trajectory  and  structure. 

The  flowfleld  over  such  general  bodies  can  be  quite  complicated.  The 
blunt  nose  generates  a  bow  shock  and  within  the  shock  layer,  depending  on 
the  body  geometry,  imbedded  shocks  and  expansion  waves  may  exist.  At  moderate 
angles  of  attack  there  could  be  a  large  region  of  crossflow  separation.  The 
flowfield  over  these  bodies  under  such  flow  conditions  can  be  numerically 
predicted  by  solving  the  steady,  three-dimensional  Parabolized  Navier-Stokes 
(PNS)  equations.  This  method  has  a  great  savings  in  computer  time  compared 
to  the  solution  of  full  Navier-Stokes  equations.  The  PNS  equations  are  para¬ 
bolic  in  the  streamwise  direction  and  elliptic  in  the  crossflow  direction 
which  makes  it  possible  for  this  method  to  solve  the  crossflow  separated  re¬ 
gion.  The  derivatives  are  written  in  the  conservative  form  which  makes  it 
possible  to  capture  the  imbedded  shocks  in  the  flowfield. 

In  recent  years  the  three-dimensional  viscous  shock-layer  approach 
(VSL3D)  *  and  the  Parabolized  Navier  Stokes  method  (PNS)^'**  have  been  applied 
to  various  problems.  The  VSL3D  equations  are  parabolic  in  both  streamwise 
and  crossflow  directions  and  are  solved  by  efficient  methods  which  require 
substantially  less  computing  time  than  the  PNS  method.  The  VSL3D*  method  can 
be  applied  to  general  geometries  to  obtain  the  flowfield  solution  over  the 
entire  body  when  the  angle  of  attack  is  not  too  high.  For  high  angles  of  at¬ 
tack,  the  PNS  methods  are  more  efficient  in  handling  crossflow  separation. 

In  this  paper,  the  general  bodies  noted  above  have  been  analyzed  by  dif¬ 
ferent  methods  and  the  results  are  compared.  The  HYTAC  code,  developed  by 
Helliwell  and  coworkers, ®  has  been  used  to  analyze  all  three  bodies.  The 
inlet  has  been  analyzed  by  HYTAC,  and  the  results  are  compared  with  those  from 
VSL3D^  .  The  results  for  the  slab  delta  wing  are  compared  with  those  from  an 
invlscid  plus  streamline  heat-transfer  code  (BLUNT  oTEIN-TRACE)®“®  and  ex¬ 
perimental  results.'® 

The  HYTAC  code  uses  a  body-normal,  shock-normal  nonorthogonal  coordinate 
system  (Figure  4).  In  this  coordinate  system  the  body  and  the  bow  shock  are 
C-C  coordinate  surfaces.  The  second  coordinate  ri  is  body-normal  and  shock- 
normal  and  always  orthogonal  to  the  C  and  ^  coordinates.  At  the  body  no¬ 
slip  boundary  conditions  are  used  and  the  enthalpy  is  specified.  At  the  shock, 
the  freestream  velocity  vector  is  transformed  into  the  computational  coordinate 
direction,  and  then  the  Rankine-Hugoriiot  jump  conditions  are  used.  The 
governing  equations  are  solved  by  implicit  differencing  in  the  n-C  plane.  A 
Gauss-Seidel  iteration  procedure  is  used  in  solving  the  equation. 


ANALYSIS 

Generation  of  a  suitable  computational  grid  is  one  of  the  important  steps 
in  the  numerical  solution  of  flowflelds.  In  the  present  study,  a  body-oriented 
nonorthogonal  coordinate  system  (C,n,C)  has  been  used.  The  first  coordinate, 

C  is  along  the  body  in  the  primary  flow  direction  (Figure  4).  The  third 
coordinate  lines  (^)  are  constructed  to  be  normal  to  the  C  coordinate  on 
the  body  surface,  and  this  facilitates  easy  calculation  of  body  forces  and 
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moments  and  easy  Introduction  of  surface  boundary  condition.  The  second  co¬ 
ordinate  (n)  is  constructed  normal  to  both  E,  and  C  -coordinate  lines  in 
the  region  between  the  body  and  the  shock  surfaces.  The  shock  surface  is 
taken  as  a  coordinate  surface  (n=l)  and  the  n-coordinate  lines  are  both 
body-normal  and  shock  normal.  A  body-oriented  cylindrical  coordinate  system 
is  used  as  the  reference  coordinate  system  for  the  interpretation  of  the  body 
geometry,  shock  shape  and  every  grid  point.  In  this  sytem,  the  ^  and  ^ 
coordinates  will  not  be  orthogonal  to  each  other  in  the  region  between  the 
body  and  shock  surfaces. 

The  steady  three-dimensional  parabolized  Navier-Stokes  equations  are 
written  in  nondimens lonal  form  in  a  general  curvilinear  coordinate  system. 

In  the  equations  for  the  stress  tensor  components,  all  the  derivatives  with 
respect  to  5  were  neglected.  The  resulting  equations  are  parabolic  in  the 
5  coordinate  direction  and  elliptic  in  the  crossflow  direction.  The  governing 
equations  are  solved  by  implicit  differencing  in  the  rj-^  plane.  The  ^  de¬ 
rivatives  are  approximated  by  a  backward  difference,  while  the  n  and  C 
derivatives  are  written  in  terms  of  an  unequally  spaced  three-point  dif¬ 
ference  scheme.  The  equations  were  linearized  by  the  Newton-Raphson  method 
and  subsequently  solved  by  Gauss-Seldel  iteration. 

At  the  surface  boundary,  no-slip  conditions  are  used  and  the  enthalpy  is 
specified.  Pressure  on  the  body  surface  is  obtained  from  the  v-momentum  e- 
quatlons.  At  the  plane  of  symmetry,  ^  =  0  and  180  degrees,  w  and  g23  are 
antisymmetric  while  all  the  other  variables  and  metrics  are  symmetric.  At 
the  shock  Rankine-Hugoniot  jump  conditions  are  used  as  the  boundary  condition. 
Since  the  shock  standoff  distance  is  also  an  unknown,  the  continuity  equa¬ 
tion  is  used  as  a  sixth  equation. 

For  a  numerical  procedure  using  a  marching  scheme,  the  construction  of  an 
accurate  initial  data  plane  is  important  in  obtaining  the  correct  solution  of 
the  whole  flowfield.  The  entire  flowfield,  coordinate  and  metrics  at  the 
initial  data  plane  have  been  generated  using  a  viscous  shock-layer  method. 

The  coordinate  and  metrics  for  the  body-normal,  shock-normal  coordinate  sys¬ 
tem  have  been  generated  from  the  two-step  body-normal  data  of  the  shock-layer 
solution  by  the  method  described  in  References  4  and  5. 

DESCRIPTION  OF  GEOMETRY  AND  FREESTREAM  CONDITIONS 

The  various  geometries  included  in  this  study  are  blunt  d  supersonic 
inlets,  an  80-deg  swept  slab  delta  wing  and  an  elliptic  body  '.his  section 

the  details  of  the  geometries  and  the  freestream  conditionc  undex  hich  they 
are  analyzed  are  explained. 

The  inlet  is  a  blunt  sphere-cone  with  an  afterbody  of  smooth  curvature. 

The  blunt  nose  has  a  radius  of  0.5  inch  and  the  cone  angle  is  12.5  degrees. 

The  afterbody  is  expressed  In  analytical  form  and  the  complete  inlet  is 
shown  in  Figure  1.  In  the  present  study,  the  external  flowfield  up  to  the 
cowl  lip  has  been  analyzed.  The  analysis  of  the  inlet  is  done  for  a  free¬ 
stream  Mach  number  of  7  at  0  and  5  degrees  angles  of  attack.  Although  the 
program  has  a  capability  to  analyze  a  laminar  flowfield,  in  this  case  turbu¬ 
lent  flow  is  assumed  and  a  two-layer  eddy-viscosity  model  is  used.  The  wall 
temperature  is  assumed  to  be  constant,  and  the  analysis  is  done  for  two 


different  temperatures.  The  details  of  the  freestream  conditions  are  given 
in  Table  1. 


Surface  heat-transfer  and  pressure  distribution  over  an  80-deg  swept 
slab  delta  wing  have  been  studied  for  a  Mach  number  of  9.6  at  zero  and  10 
degrees  angles  of  attack.  The  wing  has  cylindrical  leading  edges,  and  a 
hemispherical  nose  of  0.5-inch  radius.  The  length  of  the  body  is  about  15 
nose  radii.  The  details  of  the  geometry  are  given  in  Figure  2.  In  this  case 
the  laminar  flow  regime  has  been  considered,  and  the  surface  temperature  is 
about  540R. 

The  elliptic  body  has  a  spherical  nose  of  radius  0.3  feet  and  an  after¬ 
body  with  elliptical  cross  section.  The  change  from  spherical  cross  section 
to  the  elliptical  cross  section  is  gradual  and  is  finally  achieved  after  8  to 
10  nose  radii.  The  final  cross  section  obtained  is  a  two  by  three  ellipse. 
This  body  has  been  analyzed  for  a  freestream  Mach  number  of  20.  The  analysis 
is  done  for  laminar  flow  and  for  a  surface  temperature  of  800R,  The  details 
of  the  geometry  are  given  in  Figure  3,  and  the  freestream  conditions  are 
given  in  Table  1. 

RESULTS  AMD  DISCUSSION 


INLET  GEOMETRY 

Figures  5  and  6  show  the  variation  of  surface  pressure  and  heat  transfer 
for  the  inlet  geometry  at  zero  angle  of  attack.  The  surface  pressure  pre¬ 
dictions  agree  very  well  with  results  from  a  viscous  shock-layer  code  (VSL3D) . ^ 
The  surface  heat-transfer  predictions  by  HYTAC  is  higher  than  those  from 
VSL3D  by  about  30  percent.  Figures  7,  8  and  9  show  the  axial  velocity,  static 
pressure  and  enthalpy  profiles  for  this  case.  From  these  profiles  it  can  be 
seen  that  the  HYTAC  code  is  capable  of  capturing  the  imbedded  shocks. 

Figures  10  and  11  show  the  variation  of  surface  pressure  and  heat  trans¬ 
fer  for  angle  of  attack  5  degrees  for  various  ^  angles.  The  solution  could 
be  obtained  only  up  to  61  nose  radii,  and  it  was  found  that  at  this  point 
axial  separation  occurred.  This  can  be  verified  from  the  plot  of  axial  shear 
stress  variations  as  shown  in  Figure  12.  At  this  point  the  surface  slope  of 
tJ/e  body  is  about  20  degrees. 

SLAB  DELTA  WING 

The  computed  surface  pressure  from  HYTAC  and  STEIN  are  compared  with  ex¬ 
perimental  results  for  both  C  =  0  and  90  deg  (Figures  13  and  14).  For  both 
cases  the  results  are  in  good  agreement.  Figures  15  and  16  show  the  comparison 
of  surface  heat-transfer  results  from  TRACE  and  HYTAC  with  the  experimental 
results  for  ^  =  0  and  90  deg.  The  predictions  are  good  in  the  nose  regions 
for  both  the  cases.  For  the  downstream  region,  HYTAC  predicts  a  lower  heat 
transfer  for  both  cases.  The  predictions  by  TRACE  seem  to  be  in  error  in  the 
far  downstream  region  at  ;  =  90  deg.  This  may  be  due  to  the  fact  that  the 
computed  streamlines  diverge  from  leading  edge. 

Figures  17  and  18  show  the  surface  pressure  results  from  STEIN  and  HYTAC 
for  a  =  10  deg  for  both  i;  =  0  and  90  degree  planes.  For  ^  *  0,  the  results 
are  compared  with  the  experimental  values,  and  both  HYTAC  and  STEIN  underpredict 
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the  surface  pressure.  Experimental  results  are  not  available  for  C  =  90 
degree  plane;  however,  the  results  from  HYTAC  and  STEIN  agree  very  well. 

Figures  19  and  20  show  the  surface  heat-transfer  results  from  HYTAC  and 
TRACE  for  a  =  10  deg  for  both  ^  =  0  and  90  degree  planes.  The  results  for 
the  nose  region  agree  well  with  the  experimental  values;  however,  both  codes 
underpredict  the  surface  heat  transfer  in  the  downstream  region, 

ELLIPTIC  BODY 

Figures  21  and  22  show  the  variation  of  surface  pressure  and  heat  trans¬ 
fer  for  various  4  planes  for  the  elliptic  body.  Both  the  surface  pressure 
and  heat  transfer  are  higher  at  the  windward  plane.  Since  the  sphere-cone 
tangent  locations  depend  on  a  patching  is  needed  to  obtain  a  smooth  varia¬ 
tion.  This  procedure  resulted  in  a  slightly  concave  region  for  C  =  0  and 
30  deg  (Figure  31).  The  increase  in  surface  pressure  between  s  =  1.5  and 
2.5  is  due  to  this  concavity  in  the  geometry. 


CONCLUDING  R£>LARKS 

The  capability  of  HYTAC  in  predicting  the  three-dimensional  viscous- 
supersonic  flows  over  three  geometries  has  been  shown  by  comparison  with  ex¬ 
perimental  results  and  results  from  a  viscous  shock-layer  code  indicates  that 
HYTAC  can  accurately  solve  the  floi^jfield  over  axisynmetirc  and  nonaxisymmetric 
geometries.  Although  HYTAC  requires  more  computing  time,  it  is  capable  of 
solving  the  crossflow  separated  region.  However,  by  looking  at  the  results 
for  the  inlet  geometry,  it  can  be  seen  that  HYTAC  is  unable  to  solve  flow- 
fields  with  axial  separation.  In  this  particular  case,  even  though  the  flow 
is  attached  for  zero  angle  of  attack,  at  moderate  angles  of  attack  like  5 
degrees,  the  flow  separates  on  the  leeward  side.  In  order  to  solve  the 
separated  flowfields,  the  time  dependent  Navier-Stokes  equations  must  be  solved. 

Comparing  the  results  for  the  inlet  for  zero  angle  of  attack,  with  those 
from  VSL3D  shows  that  the  VSL3D  is  able  to  solve  the  flowfield  in  less  time. 
Although  the  surface  pressure  results  are  in  good  agreement,  the  surface  heat- 
transfer  results  from  HYTAC  is  higher  by  30  percent.  Due  to  lack  of  experi¬ 
mental  data,  the  computational  results  could  not  be  verified. 

Both  HYTAC  and  STEIN  predicted  the  surface  pressure  over  the  slab  delta 
wing  reasonably  accurately.  The  results  for  zero  angle  of  attack  agree  very 
well  with  the  experimental  results.  Both  HYTAC  and  STEIN  underpredict  the 
surface  pressure  for  10  degree  angle  of  attack.  For  the  nose  region  the  sur¬ 
face  heat-transfer  results  from  TRACE  and  HYTAC  agree  very  well  with  the  ex¬ 
perimental  results.  For  the  downstream  region,  the  prediction  by  TRACE  seems 
to  be  in  error,  and  this  is  due  to  the  fact  that  the  computed  streamlines  di¬ 
verge  from  leading  edge. 

The  results  for  the  elliptic  body  could  not  be  compared  with  other  re¬ 
sults  due  to  a  lack  of  data.  HYTAC  predicts  higher  heat  transfer  and  surface 
pressure  at  the  windward  plane.  Because  of  the  presence  of  a  concavity  in  the 
geometry  near  the  sphere-cone  tangent  point  for  ?  =  0  and  30  degrees,  the  sur¬ 
face  pressure  results  are  higher  at  these  locations. 


11-67 


In  summary,  methods  have  been  shovTn  which  can  predict  supersonic  viscous 
flowfield  over  complex  three-dimensional  bodies  at  angles  of  attack  up  to 
about  15  degrees. 
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Figure  9.  Static  Enthalpy  Profiles  for  Inlet  Geometry  at 
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Figure  18.  Surface  Pressure  Distribution  for  Slab  Delta  Wing  at 
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Figure  19.  Surface  Heat-Transfer  Distribution  for  Slab  Delta  Wing  at  a  =  10  and 
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Surface  Heat-Transfer  Distribution  for  Elliptic  Body  at 


